#! /usr/bin/perl -w use strict; my @row=(); my @seleno=(); my $ini=0; my $start=0; my $fin=0; my $length=0; my @position=(); my @ini=(); my @fin=(); my $trev=""; while () { chomp($_); @row=split(/\t/,$_); $ini=$row[8]; $fin=$row[9]; if(($row[0] eq $seleno[0]) && ($row[1] eq $seleno[1])){ push @position, $ini; push @position, $fin; push @ini, $ini; push @fin, $fin; }else{ @position=sort(@position); $ini=$position[0]; $fin=$position[-1]; @ini=sort(@ini); @fin=sort(@fin); if ($ini[0] > $fin[0]) { $trev = "-trev"; } else { $trev = ""; } $start=$ini-20000; if($start<0) {$start=1;} $length=$fin-$start+20000; print "$seleno[0] \t $seleno[1] \t $ini \t $fin \t $start \t $length \t $trev\n"; @fin=(); @ini=(); @position=(); push @position, $row[8]; push @position, $row[9]; push @ini, $row[8]; push @fin, $row[9]; @seleno=@row; } }